Note on the "Dead Zone" in Layered Accretion Models 
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V '. Abstract 

lO ! Current layered accretion models neglect the properties of the "dead zone". 

However, as argued here from simple considerations, the thickness of this zone is 
a critical quantity when the disc is in hydrostatic equilibrium. It controls not only 



the structure of the superficial, active layers, but also the mid-plane density and the 
Q\ ■ total disc mass, and should therefore be introduced in models of that kind, steady 

^ ; or not. But in the absence of intrinsic heating, the dead zone must have a tiny size 

which, given the non-stationary and turbulent character of the global flow, makes 



very likely its mixing together with the two active layers. 
Key words: accretion, accretion disks. 



Introduction 



Gammie (1996) proposed a two-phase accretion scenario for T Tauri discs to cure the 
expected inefficiency of the magneto-rotational instability (hereafter, MRI) in the cold, outer 
regions of the disc (beyond about 0.1 AU from the central proto-star) due to the low abundance 
of electrons. In this model, only the superficial layers of the outer disc are made "active" thanks 
to incoming interstellar cosmic rays (that can reactivate the MRI), leaving a non-accreting 
"dead zone" around the equatorial plane. Interestingly, the so-called layered accretion disc 
model makes two major predictions: i) an infrared excess (a common feature in the spectrum 
of T Tauri stars; e.g. Bertout 1989) caused by a positive, radial gradient of the disc accretion 
rate, and ii) some ability to develop accretion bursts (the admitted interpretation for FU Orionis 
events; e.g. Hartman, Kenyon 1996; Kley, Lin 1999) due to mass accumulation in the outer 
regions. Some aspects of layered accretion have recently been investigated : evolution of the 
solar nebula (Stepinski 1999), occurrence of eruptive events through time- dependent simulations 



(Armitage et al. 2001), linear stability properties (Reyes-Ruiz 2001), possible applications to 
accretion in active galactic nuclei (Menou, Quataert 2001) and in binaries (Menou 2002), and 
the ionization of accreted material (Fromang et al. 2002). As outlined by these authors, there 
are still many uncertainties in such a toy-model. Apart from the hypothesis (we shall not 
discuss it here) that the source of angular momentum transport is only the MRI (see Stone et 
al. 2000, and references therein), the constraint imposed on the surface density S a of the active 
layers (to be constant in space and time) appears as the strongest assumption. Misguidedly, the 
subsequent derivation of steady state profiles for the temperature, density, and disc thickness 
violates mass conservation (see also Menou 2002), or 

A7rRd t E a -d R M a ^0 (1) 

due to the non vanishing gradient of the total accretion rate, M a . Another difficulty which 
is outlined here concerns the properties of the dead zone. None of the existing models have 
yet accounted for the structure of this zone. We argue here that it plays a crucial role. In 
particular, we demonstrate from simple arguments that it is not correct to find any reliable 
solution to the problem of layered accretion (steady or not) without specifying the thickness of 
the dead zone, either in a fully ad-hoc manner (i.e. by hand), or by considering explicitly its 
structure through coupled equations, as soon as the disc is in hydrostatic equilibrium. Finally, 
it turns out that this zone should be mixed due to its small vertical extent. 

2. Status and Importance of the Dead Zone in Layered Accretion 

2.1. General Considerations 

Let us start with the basic picture of the layered accretion model as originally described 
in Gammie (1996). In this model, the disc is divided into two parts, radially: an inner region 
assumed to be fully described by the standard theory of a-discs (Shakura, Sunyaev 1973) 
and an outer, non-standard region. With respect to the central proto-star, the connection 
between the two regions is set by the temperature of gas and dust in the equatorial plane: 
below a critical value, T cr i t ~ 10 3 K, thermal ionization alone is too weak, making inefficient the 
magneto-rotational instability (MRI; see Stone et al. 2000) as the single source of angular- 
momentum transport. In the standard part of the disc, the mid-plane temperature decreases 
outwards so that the region with T < T cr i t is found at a radius larger than 1 (e.g. Frank et al. 
1992) 

i? crit -0.174 /9 q- 2 / 9 m 4 / 9 AU, (2) 

1 In this paper, all quantities and computations are scaled to a one solar mass proto-star. This choice has no 
noticeable consequences on the results exposed throughout the paper. 
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Fig. 1. Schematic view (not to scale) of the layered accretion model (left) as originally proposed by 
Gammic (1996), and expected vertical profiles (right) for the mass density, pressure, and temperature in 
the dead zone and in the active layers (see subsection 2.3). 

where a = 0.01a is the viscosity parameter (Shakura, Sunyaev 1973), is the Rosseland 
opacity of matter, 2 and m is the mass accretion rate normalized to 10~ 8 M Q yr" 1 , a typical 
value for accreting Class-II T Tauri stars (Calvet et al. 2000). In the layered disc, the MRI is 
maintained beyond -R cr i t with the help of cosmic rays, which can enhance the ionization degree 
at the disc surface, creating two superficial, accreting layers and an equatorial "dead zone" 
where gas and dust rotate without drifting inwards, as illustrated in figure 1. Assuming that 
only the MRI is present, the most critical assumption of this model is that the surface density 
of the active layers is S a m 100 g cm -2 (for each one) for R > R CT a- 

2.2. Basic Equations for the Active Layers. The Missing Closure Relation 

Let us now briefly recall the three basic equations for the disc beyond -R CT it- The first 
equation is dictated by energy balance. This equation links the effective temperature, T e , to 
the temperature 7] at the altitude z — H\ where the dead zone joins the active layer (see figure 
1). In the optically thick disc approximation, this is written as 



The second equation expresses the heating mode of the active layers. If all of the 
gravitationally energy release indeed occurs at Hi < z < Hi + h a (h a being the geometrical 
thickness of the active layer), the outgoing flux is 

2 Near 10 3 K, kr , w 3 — 5 cm 2 g~ x with a very weak dependence on the density (Pollack et al. 1994; Bell, Lin 



T ; 4 = ^ R S a T e 4 . 



(3) 



1994). 
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aT c 4 = ^ a £ a ^, (4) 

where S\ is the Keplerian rotation velocity, v a = ac a h a is the turbulent viscosity coefficient 
(Shakura, Sunyaev 1973) and c a is the sound speed. If the ionization rate has its interstellar 
medium value, cosmic rays do not contribute to the local disc heating (e.g. D'Alessio et al. 
1998). 

The third equation follows from the hydrostatic equilibrium of the whole vertical struc- 
ture. Current models neglect the dead zone at this level, but this is not correct: two equations 
must be considered with the condition of pressure balance at z — H x . Un-balancing the pressure 
would create a mean flow in the direction of the pressure gradient (i.e. vertically here). With 
the usual, vertically averaged approach and neglecting disc self-gravity (see subsection 3.3), we 
find for the active layer 

P- hT- 

fl = ^(H l + h % ) = -^-, (5) 

Pi /ia^H 

where /i a ^H is the mean mass per particle (with p a pa 2.37 for a molecular gas with cosmic 
abundances) and mn is the proton mass. A very similar equation can be derived by considering 
vertical stratification, for instance from a polytrope. This relation shows the coupling between 
the active layers and the dead zone. The point is that, physically, the thickness of the dead 
zone determines the strength of vertical gravity at z > H u and therefore the structure (i.e. 
temperature, density p a , and thickness h & = — ) of the accreting layers. Note that, with an 
appropriate formula for the Rosseland opacity, one recovers Gammie's (unsteady) solution 
from equations (3), (4), and (5) by setting H ~ (or H\ <C h a ), which ultimately means that 
there is almost no dead zone. 

2. 3. Structure of the Dead Zone 

In the spirit of the layered disc model, there is not heat generation at \z\ < H { meaning 
that the dead zone is necessarily isothermal, vertically. So, the mid-plane temperature is T m = 7] 
and the pressure and density obey a Gaussian distribution (assuming a constant chemical 
content). Actually, hydrostatic equilibrium of the dead zone yields 

/ H- \ ^ 

p m = Piexpfyj , (6) 
where A is given by [see equation (5)] 

\ 2 = ^^^2K(H i + h a ) (7) 

with /i m pa /i a (this should be correct given temperatures of interest). The surface density of 
the dead zone is then given by 

= v^ Pm Aerf (^), (8) 
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where erf denotes the Error function. It follows that a large-size dead zone (i.e. Hi ^> h & ) can 
lead to large mid-plane densities and, thus, to a large-mass disc. In other words, the dead 
zone must have a very small vertical extent compared to the active layers. This even raises the 
question of the existence of such a zone: because it is squeezed between two turbulent layers, 
it would probably mix with them and tend to disappear in a turn-over time scale. We shall 
now examine to what extent the structure of the active layers is sensitive to the presence of the 
dead zone. 

3. Essential Role of the Dead Zone 

3.1. The Strong Coupling 

We can see that equations (3), (4), and (5) involve four unknowns (T- u T e , h a , and Hi), 
leaving the system under-determined. As shown above, Hi is still unconstrained, but must be 
specified in some fashion. This is true whatever the state of the disc, steady or unsteady. In 
current models, the choice is Hi = 0, as already mentioned. An interesting question is then: 
how is the structure of the active layers modified if the parameter Hi is set to another value? In 
other words, what is the sensitivity of this special choice on the topology of layered solutions. 
To answer this question, we can simply combine equations (3), (4), and (5) to get 



which is a constant at a given radius (see footnote 2). It is therefore clear from that relation 
that Hi and h a cannot be specified independently. If h & 3> Hi, one tends to the no-dead zone 
solution and, in the reverse case, the coupling is extreme, since h\H\ ~ est. 

We can also answer the above question by considering the steady state hypothesis, which 
provides a natural closure relation for our system of four equations. In some sense, this is more 
satisfactory than imposing in equation (9) some arbitrary values for Hi, which, for instance, 
is done in Stepinski (1999). As we shall see, this hypothesis leads to an unrealistic disc, but 
this is unimportant here, since we only wish to exhibit, from a mathematical point of view, 
the sensitivity of the structure of the layered disc to Hi when a physical closure relation is 
introduced. Thus, if the flow is steady, we have (e.g. Frank et al. 1992) 



(see subsection 2.1), we must also have u a = est. This is possible with the viscosity prescription 
if cji a = est or 




(9) 




where M a is the accretion rate for the two active layers. Because S a = est in the present problem 




0.14m a 



2^-2 



AU 2 K. 



(11) 



This is the missing fourth relation. 
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3.2. Discussion of Steady State Solutions 

With the assumption of a steady state, the structure of the active layers and that of the 
dead zone can be fully determined in a self-consistent way. From equations (3), (4), and (11), 
we find that each active layer has a thickness (see footnote 1) of 

K ~ 2.6 x 1(T 2 n R 1/8 am 7/8 x 3/s AU, (12) 

where x = jjjj and the actual value of the opacity coefficient is fixed by the local temperature 
and pressure. However, the point is that h a must also satisfy equation (5) with the condition 
that Hi > 0, which is not guaranteed. Actually, equation (5) gives 

#i + /i a ~3.2x 10- 2 4 /8 "^~ 5/ V 5/8 AU, (13) 

and it follows from equations (12) and (13) that H[ vanishes for 

R m ~ 0.87K R 1/3 a" 4/3 m AU, (14) 

meaning that the steady solution exists at R > i?i n only. Generally, we have R CI it ^ R\n- Even 
i? C rit <C -Rin for usual parameter values, implying that the outer layered solution does not connect 
smoothly to the inner standard disc solution. The reason is simple: the conditions J °° pdz = S a 
and T = T crit are never satisfied simultaneously (i.e. at the same radius) in an a-disc, except if 
precisely 

rh ~ 5.4 x l(r 2 a 2 /%(T crit ) « 0.2a 2 = rh , (15) 

which is very restrictive. This smooth connection might, however, exist by relaxing the value 
of T crit or S a , for instance, or by invoking some variation of the a-parameter value with the 
radius, which is very speculative. 

Thus, for m>m , there is no steady solution between R crit and i? in . This no-solution do- 
main can be quite wide, especially for large accretion rates (as expected in very young systems) 
and/or for large a-parameter values. According to equations (12) and (13), large accretion 
rates imply large-size active layers and a small-size dead zone, whereas low accretion rates 
mean small-size active layers and a prominent, and possibly massive, dead zone. Conversely, 
for rh < rh , the total surface density at _R cr i t is less than 2E a ; the inner standard disc and the 
outer layered disc overlap between R m and -R C rit- 111 this case, an acceptable solution smoothly 
connecting the two discs should be easily found by changing the equation for the inner disc (this 
exercise is out of the scope of this paper; see subsection 3.3). But this concerns quite evolved 
(WTTS-like) systems with accretion rates lower than ~ 2 x 1O~ 9 M yr _1 (e.g. Calvet et al. 
2000), and the amount of gas available in these discs has strongly decreased (see subsection 
3.3). 

These quantitative results have been fully confirmed by numerical computation of the 
solutions 3 of equations (3), (4), (5), and (11). The results obtained for a = 0.01 and six different 

3 For this exercise, we have used realistic Rosseland opacities as compiled in Hure (2000). 
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Fig. 2. Vertically averaged solutions for the steady, layered accretion disc model (see footnote 1) computed 
for a = 0.01 and various accretion rates: rh = 0.1 (A), rh = 0.2 = rho [B; in this case, equation (15) is exactly 
fulfilled, meaning that there is a continuous connection between the inner standard disc and the layered 
disc], m = 1 (C), rh = 2.2 (D), rh — 10 (E), and rh = 13 [F; see also equation (16)]. The disc mass, Mdi sc , 
only accounts for the matter (dead zone and active layers) at R > i?; n . 

accretion rates are displayed in figure 2. There is no solution for accretion rates larger than 4 

m lim ^13« 2 , (16) 

and the sensitivity of the solution to the accretion rate and viscosity parameter value is very 
important (see the morphological change from runs C to F, whereas the accretion rate is 
increased by only a factor 6). As already announced, the steady solution appears to be non- 
physical either because it leaves a gap (the no-solution domain; see for instance run E where 
R cr[t ~ 0.6 AU and R m ~ 6 AU where T x ~ 130 K), or because it involves discs more massive 
than the central object by orders of magnitudes (see runs C, B, and A), which is not plausible. 

4 This limit is in fact compatible with the value deduced analytically by equating R m and i? ut > the marginally 
optically thin, outer disc radius where T c 7] ~ 10 K, namely i? ou t — 17m 1 / 3 AU. At 10 K, kr, — 0.02 cm 2 
g _1 (see references given in footnote 2) and the optical thickness of each active layer is ~ 2, in the Rosseland 
sense. 
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In any case, the most important point that these runs demonstrate is that the structure of 
the active layers is strongly sensitive to the structure of the dead zone, that is to Hi, which, 
indeed, can not be ignored. It is therefore expected that this sensitivity also be present in 
time-dependent models, as soon as hydrostatic equilibrium is assumed. 

3.3. Disc Self- Gravity 

The mass of gas contained in the active layers is 

M al = 4vr£ a f R R'dR' ~ 7 x 1(T 5 (JL)* M Q . (17) 

Therefore, this is a firm, lower limit for the mass of the whole disc: for instance, for 
Hi/h a = 1, the disc mass would be ~ 4 x M a i. According to Gammie's no-dead zone solution, 
hydrostatic equilibrium is expected to be modified by self-gravity as close as ~ 20 AU (the 
distance where 2^ ^ ^) from the center and beyond. In this (vertically and probably grav- 
itationally unstable) self-gravitating regime, it is possible to show again the influence of the 
dead zone on the structure of the active layers using the infinite-slab approximation (Paczyhski 
1978). Actually, equation (5) would then become 



Pi . 

— = K 

pi 



S a + is d 



(18) 



In particular, in the pure self-gravitating limit, this equation becomes 



4*Gh^ ( 1 + |f ) . 



Pi V 2£ a , 

The second term in the brackets vanishes when H r <C h a (this is no-dead zone case). 
The point is that the temperature in active layers from equation (19) is strongly sensitive to 
H^ This means that the dead zone has, again, a major influence on the structure of the active 
layers, and that the problem with massive discs invoked above remains. 

4. Concluding Remarks 

Despite its "passive" role with respect to the assumed mechanism of angular-momentum 
transport and heat generation, the dead zone is an essential component of layered accretion-disc 
models. Whatever the importance of self-gravity, this zone mechanically supports the active 
layers, and thus determines the global properties of the disc as a whole. Any layered model of 
that kind, steady or not, where the disc has reached hydrostatic equilibrium (see also Glassgold 
et al. 2000), must take the dead zone into account, and in particular its thickness, at least in 
the form of a free parameter. Probably, a model accounting for more physical mechanisms 
(and containing more degrees of freedom) could modify some issues raised here, and even 
might produce steady state solutions. However, it appears clearly that, in the absence of 
internal heating, the dead zone must have a small thickness with respect to the active layers; 
otherwise, the disc would be extremely massive. Also, without any noticeable extent, the 
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dead zone should not survive between the two turbulent layers, but mix with them. Besides, 
observations of discs around young stars indicate that they have a mass not in excess of the 
central mass (e.g. Calvet et al. 2000) as well as an outward decrease of the surface density 
(Beckwith et al. 1990; Dutrey et al. 1998), both properties of which are not compatible with 
the layered accretion model. 

It is a pleasure to thank my colleagues S. Collin, D. Gautier, F. Hersant, D. Richard, and 
J.-P. Zahn for stimulating discussions. I am grateful to the referee for important suggestions 
and comments to improve the paper. 
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